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ABSTRACT 

In a series of papers we study the stellar dynamics of the grand design barred-spiral galaxy 
NGC 1300. In the first paper of this series we estimate the gravitational potential and we give 
it in a form suitable to be used in dynamical studies. The estimation is done directly from 
near-infrared observations. Since the 3D distribution of the luminous matter is unknown, we 
construct three different general models for the potential corresponding to three different as- 
sumptions for the geometry of the system, representing limiting cases. A pure 2D disc, a 
cylindrical geometry (thick disc) and a third case, where a spherical geometry is assumed 
to apply for the major part of the bar. For the potential of the disc component on the galac- 
tic plane a Fourier decomposition method is used, that allows us to express it as a sum of 
trigonometric terms. Both even and odd components are considered, so that the estimated po- 
tential accounts also for the observed asymmetries in the morphology. For the amplitudes of 
the trigonometric terms a smoothed cubic interpolation scheme is used. The total potential in 
each model may include two additional terms (Plummer spheres) representing a central mass 
concentration and a dark halo component, respectively. In all examined models, the relative 
force perturbation points to a strongly nonlinear gravitational field, which ranges from 0.45 to 
0.8 of the axisymmetric background with the pure 2D being the most nonlinear one. The force 
perturbation in each model is found being robust to small changes of the required parameter 
values. We present the topological distributions of the stable and unstable Lagrangian points 
as a function of the pattern speed (Q p ). The topological distribution found deviates in several 
cases from the classical paradigm with two stable Lagrangian points at the sides of the bar 
and two unstable ones close to the ends of the bar. In all three models there is a range of Q p 
values, where we find multiple stationary points whose stability affects the overall dynamics 
of the system. 

Key words: Galaxies: kinematics and dynamics - Galaxies: spiral - Galaxies: structure 



1 INTRODUCTION 

Grand design galaxies are attractive objects for dynamical studies, 
because it is believed that behind their coherent global features are 
hidden the basic dynamical mechanisms that shape the morphology 
of galactic discs. Then, deviations from the grand design could be 
considered to reflect additional complications of this basic dynami- 
cal behavior. The result would be more complicated morphologies, 
like those encountered in multi-armed or flocculent spiral galax- 
ies. The motivation of this study is to investigate the dynamics 
of NGC 1300, which is a well known barred- spiral grand design 
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galaxy with a prominent bar. Its archetypical morphology promises 
a deeper understanding of the dynamics of a class of similar ob- 
jects. 

Barred-spiral systems offer many examples of grand design 
morphology, being characterized by a set of two spiral arms start- 
ing close to the ends of the bar. In some cases the spirals seem to 
emerge out of the bar as its continuation, while in others the be- 
ginning of the spirals is clearly displaced from the bar's ends (e.g. 
compare the s pirals of NGC 4535 and NGC 4548 in figure 5 in 
lGrosb0lll2OO8T) . Nevertheless the morphology of the arms in many 
of these barred-spiral s ystems is not as smooth a s that of the arms of 
normal spiral galaxies dGrosb0l & Patsislll998h . especially in near- 
infrared wavelengths. This can be characteristically seen by over- 
plotting the m = 2 co mponent of th e surfa ce brightness on images 
of galaxies, as do e.g. ISeigar et alj ([2005). One can observe that 
their arms are in general asymmetric and have gaps. It is also clear 
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Figure 1. (a) The observed K band image of NGC 1300. (b) The deprojected K band image of NGC 1300 using (PA, IA)=(87°, 35°). The galaxy was observed 
with SOFI at the 3.5 m NTT telescope at ESO La Silla. The units on the axes are in kpc, assuming a distance to the galaxy D = 19.6Mpc. The images are in 
linear intensity scale. 



that in many cases the surface brightness of the arms is larger close 
to the ends of the bar than at azimuths awa y from them. C harac- 
teristic cases are NGC 4314 and NGC 4665 fGaddotill2008h . In the 
case of NGC 4314, in near-infrar ed wavelengths, the arms hardly 
complete azimuthally a ir/2 angle lOu illen et alJl994T) . 

In the last years the dynamics of barred-spiral systems 
have been revisited by several research groups. The reason 
is the investigation of the possibility that the stars on the 
spirals are in chaotic motion (see al so articles by the same 
authors in Contopoulos & Patsis 2008 ; , and references therein). 
The motivation for this series of papers is exactly to test 
this hypothesis in one of the well known and well studied 
galaxies of this type, NGC 1300. The firs t step is to estimate 
its p otential directly from observations ids Vaucouleurs et al.l 
fl99lh . Since it is a well studied object, one can find obser- 
vations from X-rays to radio waveleng ths in order t o un- 
derstand the detailed morphology (e.g. iFabbiano et alj 1 19881 ; 
lEnglanj[l989l) . The composite color image STScI-PRC2005-01 
(http://heritage.stsci.edu/2005/01/supplemental.html ), highlights 
the characteristic asymmetry between the two sides of the galaxy. 

Images in optical wavelengths, as well as B-I color maps that 
cover an area of radius ss 6' from the centre of the galaxy, re- 
veal faint extensions o f the main bisymmetric spiral structure (e.g. 
Elme green et al.ll 996). These extensions, more clearly observed in 
bluer bands, consist of Population I objects. 

In Fig. [TJt we present the K band image of the galaxy, that 
we used in the present study. Images at near-infrared wavelengths 
depict the morphology of the old disc stellar population and are 
appropriate for dynamical studies (see Sect. [2] below). The axes 



labels are in kpc, assuming a distance to the galaxy D = 19.6 Mpc. 
At this distance l"a 0.095kpc. 

Our ultimate goal is to present dynamical mechanisms for the 
stellar component that lead to the development of the particular 
morphological features of the galaxy NGC 1300. The estimation 
of a suitable potential is the basis for all dynamical studies and 
this is done in the present paper, which is structured as follows: 
In Section 2, we present the observations we performed in order to 
obtain a reliable estimation of the potential of the galaxy. In Section 
3 we describe the image processing techniques that allowed us to 
create images that will be used to compare our models with. The 
estimation of the potential is given in Section 4. The forcing in 
the models is presented in Section 5, while the properties of the 
corresponding effective potentials in Section 6. Finally we discuss 
our conclusions in Section 7. 



2 OBSERVATIONS 

The distribution of luminous matter in a galaxy is indicated by its 
surface brightness map. The translation from luminosity to mass, 
through a mass-to-light ratio (M/L), may vary significant depend- 
ing on both the spectral band used and the underlying stellar pop- 
ulation. The near-infrared K band at 2.2 fi is close to the emission 
peak of the old stellar disc and bulge populations whi ch represent 
the m ajor visual mass constituency of a spiral galaxy dRix & Riekd 
Il993h . Although very young stars and red super giants also con- 
tribute to the K band flux, it remains a good indicator of the mass 
distribution of luminous matter outside star forming regions. 
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The K-map of NGC 1300 was observed on 2002-09-01 with 
SOFI at the 3.5 m NTT telescope, ESO La Silla. The total expo- 
sure time on target was 10 min in the K s filter, at 2.162^, with 
equal time for sky fields. A jitter pattern with 10" shifts was used 
to allow removal of bad pixels through stacking, while sky frames 
were interleaved with offsets around 10' from the galaxy. The SOFI 
instrument had a Rowckwell Hg:Cd:Te 1024x1024 detector with 
0.29" pixels on the sky. The field covered the main spiral struc- 
ture of NGC 1300 but did not provide an accurate estimate of 
the sky background level. The final K band frame had a seeing of 
0.8" and reached a surface brightness of K = 20.8 mag arcsec~ 2 at a 
signal-to-noise (S/N) level of 3. Thus, the quality of our dedicated 
SOFI/NTT observations can be considered as an additional advan- 
tage for potential calculations, besides the known advantages of 
the K-band for studies related with the mass distribution in galax- 
ies. For comparison, the 2MASS survey offers a seeing of 2-3"and 
the corresponding images reach 20 mag arcsec~ 2 at a S/N level 
of l<x (see pages under http://www.ipac.caltech.edu/2mass/!. The 
reduction of the SO FI data followed the procedure described by 
lGrosb0letal.1 ( l2OO4h . 

The sky proje ction parameters were adopted from 
iLindblad et alj d 19971) who estimated the Position (PA) and 
Inclination Angle (IA) by fitting tilted ring models to their HI 
velocity data. The parameters values they found were (PA, IA) 
= (87°, 35°). These values are different from those one finds 
by fitting an exponential disc to the regions outsid e the main 
bar in the K band image, as in lGrosb0l et all {2004), which are 
(P A, IA)=(106.6°,42.2 ° ). Th is values are close to those found 
bv lLaurikainen & Said Hool) . However, since for barred spirals 
like NGC 1300, any determination of the sky projection based 
on surface photometry may be strongly biased due to the strong 
m = 2 mode and the open spiral stru cture, we preferred to use 
for our study the values proposed by ILindblad et all dl997l) . A 
disadvantage of the use of the near-infrared image is that one 
can hardly detect the axisymmetric disc in the outer regions. 
We have to note though, that for the response models and the 
orbital analysis, we investigated also the effect of the deprojection 
parameters in the potential, and the results are presented in the 
forthcoming papers of this series. 

The deprojected with parameters (PA, IA) = (87°, 35°) image 
of NGC 1300, is shown in Fig.[Tp and depicts the morphology we 
will try to model. The size of this image is 531x531 since it has 
been re-binned to 0.5"pixel size. The standard background subtrac- 
tion for the K-band frame could be slightly in error since it was 
based on separate sky exposures and not on the frame itself which 
was almost fully occupied by the galaxy. 



3 DATA REDUCTION 

We adopted as distance of the galaxy the value D = 19.6 Mpc, 
according to the cosmology-corrected luminosity distance given in 
the NASA/IPAC Extragalactic Database (NED). This value may 
differ from the distance used in other papers on the dynamics of 
NGC 1300. In any case, the exact distance of a galaxy does not 
influence essentially the dy namical phenomena that shape its mor- 
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Before the K-band frame can be used to bootstrap the esti- 
mate of the potential, two possible corrections should be consid- 
ered namely the population effect on the (M/L) ratio and the exact 
background level. Although the major fraction of fl ux in the K- 
band originates from the old stellar, disc population dRix & Riekd 



Il993h . young stars in the arm regions will bias the effective M/L. 
The young stellar population may be present both as a diffuse, un- 
resolved, component and as resolved young ste l lar complexes seen 
in the arms as bright knots jPatsis et al]|200ll ; lGrosb0l & Dottoril 
2008). In order to reduce the M/L bias from young stars, we re- 
moved such bright knots by a digital filter. A correction that has to 
be done concerns the brightness of the outer parts of the galactic 
disc. This area is biased towards higher values due to the presence 
of young stellar c lusters and stellar complexes, still present even 
in K band images jPatsis et al.l2QOll ; lGrosb0l & Dottorill2008l) . For 
this purpose, in Fig. [Tp we firstly removed the bright spots in the 
region with a radius R > 6.5kpc. Bright spots, assumed to be young 
clusters, do not cover areas larger than 10 x 10 pixels. We substi- 
tuted the intensities of all extreme high (but also a few low) valued 
pixels with the median value of the corresponding 961 (31 x 31) 
pixels area centred at the very pixel under consideration. In Fig.[2ji 
only the identified as high valued pixels are shown, while in Fig.[2p 
we see the rectified image of the galaxy after the removal of the 
bright spots. After these reductions the galaxy image seems quite 
smooth in the area of high valued pixels (Fig. |2p) but still presents 
significant noise in the area of low valued pixels. In Fig. [2J; we 
see how noisy is the contour corresponding to the brightness value 
20.8mag arcsec~ 2 . In order to smooth out the image, we apply a 
low-pass smoothing filter which affects only the low valued pixels 
(> 19mag arcsec~ 2 ) getting a less noised image (Fig.[2ji). 

In Fig. |2jl we indicate four radial directions A, A', B and B', 
centred at the axes origin, along which we plot the surface bright- 
ness after applying the smoothing filter. We have chosen the 
directions along the minor axis (A, A') which reveal better the ex- 
ponential disc. The directions (B, B') are along the diagonal of 
the image nearest to the minor axis and provide data for longer 
distance. The surface brightness along these directions is given in 
Fig.[3ji. It is obvious that no exponential profile is observed at the 
outer parts of the disc (beyond the region with the bar and spiral 
structures) as expected. This implies that the adopted zero value 
level for the background is not correct. Therefore, we decrease the 
zero value level by the smallest quantity so that the surface bright- 
ness profiles become exponential, i.e. 



Z(R) oc e 



■R/hg 



(1) 



as shown in Fig.[3p (actually we increase the brightness values of 
all pixels by the same quantity). We consider the region with the 
spiral arms part of the exponential disc. We note that the applied 
correction is of the order of the background uncertainty. The cor- 
responding scale length we find in that case is h R ^ 10 kpc. By 
increasing the brightness by a smaller value than the one we finally 
applied, we do not get a well defined exponential profil e. The value 
of 10 k pc we find is in agreement with the one given bv lPrieto et al.l 
(2001) corresp onding to observatio n s of N GC 1 300 in the I band. 
We note that lLaurikainen & Said d2002h and lLaurikaine n et al.l 
d2004l) find values close to 7 kpc in their 2MASS and H-band im- 
ages respectively. The effect of such differences on the nonlinearity 
of models for the potential is investigated below (Sect. 5). 



4 THE ESTIMATION OF THE POTENTIAL 

After all the above amendments we result in an estimation for the 
surface density distribution of the luminous component in the re- 
duced image (Fig. |2p), under the assumption of a constant mass- 
to-light ratio. In our approach we reconstruct the galactic potential 
assuming the maximum contribution of the luminous component 
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Figure 2. (a) The bright spots identified as young clusters. These regions correspond to extreme high valued pixels, (b) The K band image of NGC 1300 after 
reducing the intensities of the extreme high valued pixels, This image is smoother than the one in Fig.[T]) and lacks of bright areas corresponding to young 
clusters, (c) The contour line for 20.8mag arcsec~ 2 for the image panel (b). It is obvious that the low brightness data are noisy, (d) The same contour line, as 
in (c), after the application of a low-pass smoothing filter to low brightness pixels (> 19mag arcsec~ 2 ). The arrows A, A', B, B' denote the directions along 
which we plot the surface brightness in Fig. [5] 



Table 1. All the parameter values used in j7}, {U and i 1 6i for each model (see text). These values correspond to the rotation curves of the models in Figs.|4]El 
and[5] The last two columns read the total mass inside the (spherical) radius r = 15kpc for the luminous (LC) and dark (DC) component respectively. 



models 


M C mc(M ) 


"CMC (kpc) 


M DH (M Q ) 


a DH (kpc) 


io(M kpc 3 ) 


umhp (kpc) 


LC (M Q ) 


DC (M ) 


Model A 


7.0 X 10 s 


0.1 


8.0 x 10" 


30 






5.0 X 10 10 


7.6 x 10 10 


Model B 


2.3 X 10 9 


0.1 


1.6 x 10 12 


40 






6.5 x 10 10 


7.5 x 10 10 


Model C 


9.4 X 10 8 


0.07 


1.6 x 10 12 


40 


6.12 x 10 10 


0.19 


7.1 x 10 10 


7.5 x 10 10 



to the observed rotation curve as given in iLindblad et ail Jl997t) . 
For the cases that parts of the rotation curve will not be able to be 
entirely assigned to the luminous component we switch on two ad- 
ditional components in our modelling. They represent a dark halo 
and/or a central mass concentration. Since our goal is the study of 
the dynamics of the stars, treated as test particles moving in the es- 



timated gravitational field, the total potential must be compatible 
with the rotation curve of the galaxy, in every model we study. 

In order to proceed with the potential estimation of the lumi- 
nous component we need to do an assumption about the distribution 
of the matter in the third dimension. Below, we study three differ- 
ent cases. The basic difference among these three general models 
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Figure 3. (a) The surface brightness along the directions indicated in Fig. |2ji. We observe that no exponential law can be identified in the outer parts, which 
means that the considered zero level is not correct, (b) The surface brightness along the same directions after the minimum correction of the zero level so that 
an exponential decrement is obtained in the outer parts. The scale length of the resulting exponential law is 10 kpc. 



is their geometry. We want to compare the marginal thin geometry 
of the 2D case (model A) with a thick disc case (model B) in order 
to see if the comparison between galaxy and models improves by 
varying our parameters towards the one or the other direction in the 
parameter space. Model C changes essentially the geometry of the 
bar from cylindrical to a combination of spherical with cylindrical 
and should be again considered as a limiting case. More specifically 
we consider: 

Model A. The '2D' or 'degenerate' case in which all matter is 
considered lying on the (x, y) galactic plane. 

Model B. The 'disc' case in which we assume a cylindrical ge- 
ometry and that all matter is distributed in a 3D disc with scale 
height h z . Then the vertical dependence of the density p reads 
jvan der Kruii 19881) 



1 



p<xp z (z) = rr-sech 
2h- 



2/i, 



(2) 



where h z is the vertical scale height and n is an index (see Section 
4.2, for its meaning). 

Model C. The 'spheroidal' case in which we consider that the 
total observed brightness is the combination of two major com- 
ponents: i) a spheroidal component that accounts for the ma- 
jor part of the bar and ii) a 3D exponential disc with scale 
height h z as in mod el B. The spheroidal component includes the 
bulge (r < 1.5kp c 



:1 r>. the spheroidal component includes the 
Prieto etal.|[200ll : lLaurikainen & Saldl2002l : 



lLaurikainen et alj|2004l) and the axisymmetric contribution of the 
central part of the bar. 

In Cartesian coordinates for a given spatial density p(x,y, z) 
the gravitational potential 0(x, y, z) is 



J V(*-x') 2 + (y-y') 2 + (z-z' 



) 2 + (y-y') 2 + {z-z') 2 



(3) 



where G is the gravitational constant. This means that 9(x, y, z) 

is the convolution of the functions p(x, y, z) and g(x, y,z) = — = 

r 



Thus the potential is written as O = p ® g, where < 



V* 2 +r + z 2 

denotes the convolution operator. 



4.1 Model A 



In model A, where all mass is assumed on the (x,y) plane, the den- 
sity can be written as p(x,y, z) = X(x,y)<5(z), where S is the surface 
density and 6 is the delta function. Substituting the above density 
expression into ([3} and integrating it over z we get 



®(jc,y,0) : 



; Jz(x',y' 



)g(x - x',y - y, 0)dx'dy' 



(4) 



which is the convolution of the functions E(x, y) and g(x,y, 0). 
Knowing the surface density 'L{x i ,y l ) at each grid point (i, f) we 
can calculate the potential at the same grid points by using the con- 
volution theorem 



<d = r- [ [T(2) rig)] , 



(5) 



where T and T~ 1 , denote the Fourier transform and inverse Fourier 
transform respectively. We use an extended square area 16 times 
larger than the original area of the image of the galaxy in order 
to avoid the contribution from periodic reflections of the signal 
from the galaxy's pixels. Such an extended area will allow us to 
consider, in dynamical studies, stars that travel to large distances 
away from the centre of the galaxy. The surface density in the ex- 
tended area is considered to be zero. We use the 2D fast Fourier 
transform (FFT) technique for the corresponding computations. 
Note, that in this calculation the adopted g(0, 0, 0) value is equal 
to (bin size) -1 ^21 kpc~' since the softening length is of the order 
of the bin s ize. A detailed de scription of the whole procedure can 
be found in lHohl & Hocknevl fl969t) . 

Given the potential values at the grid points one can calculate 
the potential everywhere in between these points using a 2D inter- 
polation scheme. However, such a scheme, besides its complexity, 
provides the potential taking into account non-smoot h, local den- 
sity disturbances. For this reason we proceed as in lOuillen et al.l 
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Figure 4. The rotation curve of model A together with observational data by 
Lindbla d et all <1997t) . Gray dots correspond to optical measurements from 
various slits aligned close to the major axis of the galaxy. The gray curve 
with the enor bars conespond the HI rotation curve given in Lindb lad et alj 
il997t) . The dotted line represents the rotation curve of the luminous com- 
ponent only along the major axis of the bar in model A. The dashed curve is 
obtained by considering also an additional central mass (CMC) and a dark 
halo (DH) component. The importance of these terms are determined by the 
parameters given in Table 1 . 



and we decompose the potential into its azimuthal Fourier 
components at every radius R. Under this consideration the poten- 
tial can be written in polar coordinates (R, 8) as 

3>(R, 8) = <t> (R) + £ [<b kc (R) cos(W) + sin(fee)] . (6) 

For kmax — * °° the form ® returns exactly the same values with 
those we have from the 2D interpolation scheme mentioned above. 
We cut off the high harmonics and consider all terms (even and odd) 
up to k max = 6. Like this we are able to reproduce all the important 
morphological featu res of the galaxy (its asymmetries included). 
Ouille n et all d 19941) calculated the amplitudes O fc (tf), <3>t<(R) for 
many rings R f and fitted them by 8th order polynomials. Despite 
the advantages of simplicity and analyticity of the expressions for 
each function that the use of polynomials secures, it is not capa- 
ble of following the actual variations well. For this reason in the 
present study we apply the interpolation scheme with cubic splines 
to the (R h <b kc (Rj)), (R h O tt (i?,)) data sets. This scheme guarantees 
the continuity of the potential, of the forces and of derivatives of 
the forces, which are used in the variational equations providing 
indicators that will measure the chaoticity of the orbits (e.g. Lya- 
punov exponents). This is particular important for stellar dynami- 
cal investigations, which is our motivation for the estimation of the 
NGC 1300 potential. Note that we set O fc (0) = O t [ (0) = in order 
to secure the continuity required at R — > 0. Table 2 in Appendi)Q 
presents all the amplitude values, at 101 successive radii Rj, used 
for the derivation of the interpolating splines. 

Figure[4]shows the HI rotation curve (gray solid line with er- 
ror bars together with optical velocity measurements (gray dots) 
from variou s slits aligned close t o the major axis of the galaxy, as 
taken from Lindblad et al. Ul997l) . In this figure we have also plot- 
ted a rotation curve along the major axis of the bar, corresponding 
to the potential {6]l of the luminous matter only (dotted line). In this 
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case it has been considered being 5.0 x 10 10 M o (see Tablefj}. This 
curve fails reproducing the flat part of the rotation curve at large 
distances. This implies that some additional mass component exists 
besides the luminous component. The most reasonable candidate 
to be invoked is a dark halo (DH) term. As already mentioned our 
general model foresees the switching on of two additional compo- 
nents for a central mass concentration (CMC) and a DH. We model 
both of them with a Plummer sphere since it is simple, effective and 
widely used in the literature. The importance of this terms are con- 
trolled by the parameters of the Plummer spheres. Thus the total 
potential <t> T (R, 8) can then be written 

4> T (R, 8) = <f LM (R, 8) + <t> CMC (R) + <S>dh(K) (7) 

where O tM is the luminous matter component given by l|6} and 

GMcmc 

®cmc(R) = DH - (8) 

DH IR 1 + a 2 

/n t- u CMC 

\ DH 

is the potential for the Plummer spheres, with a being a constant 
parameter. The adopted parameters are given in Table[T](Model A). 
In Fig. 5 we see that the DH term is quite important and helps the 
rotation curve of the model reaching the level of the HI curve of 
Lindblad et al. (1997). On the other hand the CMC term is less im- 
portant increasing the central peak by ss 12km/sec which eventually 
reaches the value as 208km/sec. Note that the central peak value is 
not well defined by the observational data and therefore the adopted 
parameter values for the CMC term provide only a reasonable area 
in the parameter space. 

4.2 Model B 

In case B we assume that the density can be written as 

p(x,y,z) = V{x,y) Pl {z). (9) 
Ivan der Krui I dl988l) proposed a family of models regarding the 



vertical dependence of the density which reads 



p z (z) = Asech ^— j with n e [1, oo) 
where the normalization constant A is given as 



1 + ;) 



>that£"p z (z)rfz 

r(w) 

Varying n from n = 1 to n 



2/*,Vir(i + i)' 

1. r is the Gamma function defined by 



l e-'dt. 



(10) 



(11) 



(12) 



co we get p z (z) varying from 



p z (z) oc sech" 



(isothermal disc) to 



p,(z) oc exp(- — 



2h 



\z\ 



(exponential disc). Ide Griisl \ 19981) and lKregel et ail ( 2002) provide 
estimations for the mean values of the parameter n and the ratio 
h R /h z in samples of galaxies, which ar e n x 3.7 ± 0.6 and h R /h z « 
7.3 ± 2.2, respectively (for a review see lvan der Kruitll2002h . These 
studies have been based on observations in I and K bands. 

We adopted the values n = 4, h R /h z = 7, which give 
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Figure 5. The rotation curve along the major axis of the bar for model B. 
The drawn curves are as in Fig. [4] We observe that the luminous component 
alone (dotted curve) reaches lower velocity values at the central and outer 
parts of the galaxy with respect to the velocity measurements. The added 
CMC component is a 3 times more massive than in case A. The rotation 
curve corresponding to the total potential {7) (dashed curve) with the pa- 
rameter values shown in Table[T]levels off at about the same u max as the HI 
rotation curve. Note that the central peak is sharper in this case. 



h z ss 1 .43 kpc and following lOuillenet all Jl 9941) . we substitute the 
density expression ([9} into and by integrating over z we obtain 



®(x,y, 0) 



E(x',/)g B (x - x ,y - y')dx'dy' , 



where the function g B is given by 



gis(x,y) ■■ 



p-Xz')dz' 



(13) 



(14) 



y + Z '- 



Thus, the potential O on the galactic plane is the convolution of the 
functions E and g B . We can again calculate it applying the convo- 
lution theorem O = IF -1 [F(X) T{gB)\, since we know the surface 
density £(x,-, yf) at every gri d point (i, /'). 

We note that lLaurikainen & Said d2002h and 
lLaurikainen et al.l d2004h introduced an improved version of 
the Quillen et al (1994) method, where the image pixels are first 
smoothed by calculating the azimuthal Fourier decompositions of 
the surface density in different radial zones in polar coordinates. 
This approach has some advantages, like a possibility to use more 
noisy images and to apply radially non-constant vertical models. 
However, for high quality images and for strongly barred galaxies 
(like our data for NGC 1300) the differences in the resulting 
potential should be negligible regardless of the integration method. 

Working in the same way as in model A we obtain expres- 
sions mathematically similar to 01 for the potential correspond- 
ing to the luminous component and to for the total potential. 
Table 3 in Appendix contains all the values of the amplitudes for 
<S> (Ri), (frkciRj), 9 ks (Rj) calculated for 101 rings at Rj. The specific 
values for ® > ®fa correspond to 6.5 x 10 10 M o for the mass of 
the luminous matter. 

Finally, Fig. [5] presents the rotation cu rve corresponding to 
model B together with the observational data (Lindbla d et al J 19971) 
as we did in Fig.|4]for the 2D case. The dotted line gives the rota- 
tion curve along the major axis of the bar when we take into ac- 
count only the luminous component. We observe that this rotation 
curve reaches lower values than the observational one also at the 
central parts. For this reason we add a CMC term which is = 3 




lit 12 14 



R(kpc) 



Figure 6. The surface brightness of the MHP term in model C. The solid 
line represents the mean (over azimuth) surface brightness after the sub- 
traction of the exponential disc surface brightness (see Fig. [3p) from the 
total one. This part is fitted by the MHP with the parameter values shown in 
Table[JJ(dashed curve). 



times more massive than the one we added in model A, i.e. now 
it is 2.3 x 10 9 M o . We add also a DH term in order to get a bet- 
ter description at the outer parts of the disc. The dashed line gives 
the rotation curve of the total potential with the corresponding 
parameter values given in Table Q] 



4.3 Model C 

In this case we consider that the main part of the bar can be de- 
scribed with a spherical component. This part could be vaguely de- 
scribed as "the bar without its ansae" (Fig. []])). For this purpose we 
use th e modified Hubble profile (MHP) (see iBinnev & Tremaind 
I2008L p. 68) 



Pmhp - Jo 



1 + 



-3/2 



(15) 



The parameter a MHP is the scale length and the parameter jo deter- 
mines the mass inside a specific radius (e.g. r < a MHP ). Note that 
the mass inside a radius, in MHP model, diverges logarithmically 
at large radii and therefore it is meaningless to refer to total mass. 
The corresponding potential via Poisson equation is 



#1 



GM,,(r) AnGj Q a 2 Mi 



\ V "MHP I 



(16) 



where 



M h (r) = AnalfupjoX 
r 

In 



OMHP 



1 + 



<*MHP 



O-MHP 



1 + 



®MHP 



-1/2 



(17) 



By integrating d 1 5 b over z we take the expression for the sur- 
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Figure 7. Rotation curve together with observed velocities for model C. 
Drawn curves are as in Figs. [4] and \5\ The rotation curve of the luminous 
matter along the bars' major axis is given with the dotted curve. The corre- 
sponding curve of the total potential including CMC and DH components, 
with the parameter values given in Table [TJ is represented by the dashed 
curve. 



face density for MHP, which is 



o(R). 



2 , p2 



(18) 



In order to find the values of the best fitting parameters of the MHP 
model we subtract the exponential disc l.j oc e~ s /*« with h R = lO kpc 
(see Fig.(3j)) from the observational surface density Y. obs so that the 
remaining luminosity at the central region (R < 4 kpc) is mainly 
due to the assumed spherical component. In Fig. [6] we plot the az- 
imuthally mean surface brightness of the central area of the galaxy 
(solid line) after the subtraction of the exponential disc. We also 
give the curve corresponding to J 1 8t (dashed curve) with the pa- 
rameter values shown in TableQ] We observe that the MHP gives a 
reasonable fitting of the inner 3.5-4 kpc. 

The total potential corresponding to the luminous mass is the 
sum of the MHP component l |16t and the component corresponding 
to the remaining mass. The surface density of the latter component 
So is derived by subtracting the surface density TLmhp from the sur- 
face density E o6s we derive from the observations. Thus, the surface 
density (of the remaining disc) Ed(x,, yf) at each pixel (i, j) is given 
by 



^D(Xi,yj) = Z, obs (Xi,yj)- < Z, MHP (x h yj) >, 



(19) 



where < l. M Hp(xj,yj) > is the mean surface density ~L M hp inside 
the area of the squared pixel centered at (x h yj). We consider this 
component as a 3D exponential disc with vertical distribution given 
by d 1 Ob - In order to obtain an expression similar to ((SJ for the disc 
in model C we work in the same way as in model B. Table 4 in Ap- 
pendix gives the values of the amplitudes of ® WX ^fa-C^i), W) 
of the potential corresponding to the disc component. 

Figure [7] shows the d ata of the observational rotational veloc- 
ities dLindbladet ai]|l997h as in the two previous models together 
with the rotation curve of the luminous component (both MHP and 
disc) along the major axis of the bar represented by the dotted 
curve. Also in this case the activation of the two additional poten- 
tial terms, especially the DH one, allows a better relation between 
the rotation curve of the model along its bars major axis and the 
inclination of the observed HI curve at large distances. The values 
of their parameters for model C are also shown in TableQ] 



5 THE NONLINEARITY OF THE MODELS 

As a measure of the nonlinearity of the models we consider the 
quantity 



Q 



|AF| |F to tal - F ax i s | 



|F a xisl 



|F ax isl 



dR 



dR 



, 1 d<b Ma X 



(20) 



dR 



which is the ratio of the modulus of the vectorial difference of the 
total force (axisymmetric and non-axisymmetric) from the axisym- 
metric force, over the modulus of the axisymmetric force. Thus, 
Q measures the non-axisymmetric force perturbation, normalized 
over the axisymmetric one. This parameter describes well the non- 
linear effect of the non-axisymmetric components but it is unable 
to determine by itself the global dynamics in a rotating system (see 
below the discussion about the importance of Q. p ). 

In the past, the strength of bars and spirals i n disc galaxies 



has b e en quantified by m eans of similar indices (e.g.y 
l200lUVorobvovl l2006h . The index introduced by 
200 lh has been w idely used in large samples of galaxies 



Buta & Block! 



Buta & Block! 



Laurikainen & Said 120021 : iBlock et alJ 120041 : lLaurikainen et al.l 
l2004h . This index is similar to the quantity Q defined in Eq. ( 1201) . 
The variation of Q in our models is presented in Fig. [8] 

At the left-hand column the panels show in color diagrams 
the strength of the non-axisymmetric force, measured by Q, for the 
three models. Each model is indicated at the right edge of the fig- 
ure. At the right-hand column we plot the maximum Q as a func- 
tion of R (black solid curves). We observe that the overall maxi- 
mum perturbation strength reaches higher values (« 0.8) in model 
A (pure 2D), while in the other two models it is lower reaching 
x 0.45 and « 0.5 for models B and C respectively. The maxi- 
mum Q values of models B and C are close to those presented by 
Laurikaine n~& Said |2002), while the pure 2D case has a significant 
larger overall maximum Q value. In each panel of the right-hand 
column we have plotted additional maximum Q curves correspond- 
ing to different parameter values, which are indicated in the figure. 
The green curves give the maximum Q profiles without the DH 
component. We observe that the maximum Q values, with the DH 
component included, are approximately 0.8-0.85 of the correspond- 
ing values neglecting the D H component. Suc h values are within 
the range of values given bv lButaetai1 ( l2004l) . although towards 
the lower limit. We studied also the robustness of the maximum Q 
profile to the adopted orientation parameters. We considered two 
additional orientation parameter sets i.e. (PA, IA) = (106.6°, 42.2°) 
(dashed magenta curve, Fig.[8p) and (PA, IA) = (67°, 28°) (dashed 
blue curve, Fig.(8j)). The first set of the o rientation parameters is the 
one obtained from near-infrared data bv lGrosb0l etaTJ J2004h . The 
second set has angles deviating from the values we finally used as 
much as the first one, but to the opposite direction. We observe, that 
the maximum Q value, in these cases, is affected by less than 10%. 
We also checked the dependence of the maximum Q value on the 
changes of the zero-value level, which yield to disc scale lengths 
of 7 kpc (dashed red curve, Fig.[8p) and 13 kpc (dashed-dotted or- 
ange curve, Fig.[8p), respectively. Note, that the exponential discs 
with these scale lengths do not fit the corresponding data as well as 
in the original h R = \Q kpc case. We observe, that the maximum Q 
value, in these cases, is affected by less than 5%. 

In Fig. [8]i we see that the maximum Q value, in model B, 
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Figure 8. The Q variation in our three models for several sets of the potential parameters. The panels of the left-hand column give Q for models "A" (a), 
"B"(c) and "C"(e). The panels of the right-hand column give the maximum Q as a function of R (black solid curves). Additional curves refer to the radial 
variation of maximum Q in models with different parameters than those used in the basic set-up of the three general models. For model "A" in panel (b) the 
dotted green curve refers to a model without dark halo, the magenta and blue dashed curves to different orientation parameters, while the dashed pink and dark 
yellow ones to models with disc scale lengths Iir = 1 kpc and Iir = 13 kpc respectively. Arrows help to better trace each curve. For model "B" in panel (d), 
additional curves are given for a "model B" without dark halo component (dotted green curve) and models with hnjh- =5,9. Finally for model "C", panel (f), 
an additional curve gives the radial variation of Q again in a model that with respect to "C" does not have any dark halo component. 
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Figure 9. The Q value maps for model A on the galactic plane for k lnax = 6, k max = 12 (first and second row respectively) with and without odd terms 
(left-hand column and right hand column respectively). A combination that can be used for reliable dynamical models is with k max = 6, including all terms 
(panel (a)). 



varies approximately by 10% when the ratio h R /h z varies ltr from 
its mean value (see dashed-dotted orange and dotted red curves 
in Fig. [8}l). We also certified, that the Q profiles almost coincide 
within lcr variation of the exponent n of Eq. d 10b from its mean 
value. 

Special attention has been given to the number and kind (even 
or odd) of the Fourier terms to be included in the potential. The 
morphological asymmetries observed in the galaxy underline the 
importance of the inclusion of the odd terms for dynamical stud- 
ies. The number and the kind of the Fourier terms may affect the 
maximum Q value as well as the 2D distribution of the Q value on 
the galactic plane. In order to study these effects we constructed Q 
value maps on the galactic plane for different combinations of the 
Fourier terms (Fig. [9j. We present the results for model A since it 
is the model with the highest non-axisymmetric amplitudes in the 
Fourier decomposition. The upper two panels (a,b) of Fig.[9]include 
terms up to k max = 6 in Eq. (SJ, while the two panels at the bottom 
(c,d) up to k max = 12. The left panels (a,c) include all k terms, 
while the right ones (b,d) only the even terms, as indicated at the 
top of each of them. By including all terms, the maximum Q value 



for the case we consider up to k max = 6 is 0.96 of the maximum 
value corresponding to k mas = 12. By including higher order terms 
(kmax > 12), the maximum Q value does not get further refinement. 
On the other hand, by neglecting the odd terms, the Q map clearly 
fails to reproduce all the asymmetries and the maximum Q value 
becomes at least 10% lower than that with the odd terms even for 
the k max =12 case. More specifically, the maximum Q values cor- 
responding to panels (a), (b), (c) of Figure|9]are 0.96, 0.89, 0.84 of 
the maximum Q value we get in Fig. [5J:. We can observe, that at 
the sides of the bar the prominent green "tails", which essentially 
describe the force field at the spirals' region, are conspicuous only 
in panels (a) and (c). Also at radii around 8 and 9 kpc, at the ends of 
the bar, only panels (a) and (c) give a rather constant green shade, 
while at the same radii green color fades out in panel (b), but also in 
(d). Evident is that in terms of nonlinearity of the potential the odd 
terms are essential, while including k > 6 terms does not improve 
the models at a level worth the additional computing time needed 
for the calculations. Thus, a reliable combination that can be used 
in dynamical studies is k max = 6 with all even and odd terms in- 
cluded (Fig.|5Ji). 
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Figure 10. The loci of the Lagrangian points (L[ , L2 in the left-hand column and L4 , L5 in the right-hand column) for all the models as a function of the fi p . 
The red (blue) points correspond to unstable (stable) Lagrangian points. In each panel we have projected the Lagrangian points on the galactic plane. We 
observe that in all models there are ranges of Q„ values for which we have multiple Lagrangian points. This property is more evident in model A, where we 
observe longer ranges of Sl p values corresponding to multiple Lagrangian points. Note that the separated group of points which are located at high x values 
(seen in all the left-hand column panels) correspond to multiple L[ Lagrangian points ordered near-parallel to the y axes and they refer to cases with low fi p 
values. 
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Figure 11. Some characteristic examples of multiple Lagrangian points. We plot the contour lines of the effective potential overlapped with the image of the 
galaxy for the cases and the Q„ values indicated in the figure. Big yellow (green) dots correspond to stable (unstable) Lagrangian points. In (a) we see 3 
unstable and 2 stable L\ points. In (b) we have 2 stable and 1 unstable L5 points. In (c) finally we see 2 unstable and 1 stable L2 points. Note that the considered 
sense of rotation is clockwise. 



In all examined cases the forcing of the models suggests 
that our systems are strongly nonlinear, since the perturbing term 
is at least higher than 40% of the axisymmetric background. 
This is a usual situation in m odels of strong barred galaxies 
dKaufmann & Con topoulos 1996). 

We note that potential evaluatio n s for NGC 13 hav e 
been previously done by lEnglanoj (fl989l) . iHunterl Jl990t) . 
iLindblad & Krister] JT 996). Ague rri et al J d200 it) and more recently 
by Laurikainen & Salo ( 2002) using a 2MASS K-band image and 
bvlLaurikainen et al. ( 2004) by means of OSUBGSG H-band data. 
These potentials have some features similar with our general model 
(for a discussion see Sect. 7). In our case, emphasis is given in con- 
structing a potential function suitable for dynamical studies. We 
kept fixed the quantities that could be derived directly from our 
near-infrared K-band data, while the parameters that could be only 
confined within some range have been treated as free parameters. 
However, the most important parameter for the dynami cs of a rotat- 
ing d isc galaxy is the pattern speed (for a review see lConto poulos 
2004), which is an unknown quantity. The most simple assumption 
is that the system we are studying has a unique pattern speed. We 
take it as a basic hypothesis, the validity of which will be evaluated 
from the feedback we will get from response models (to be pre- 
sented in subsequent papers). At any rate, it is the effective potential 
that determines the force field in our calculations. The dynamics of 
two systems with the same potential O but with consi derably dif- 
feren t pattern speeds are considerably different (see e.g. lPatsis et al.l 
Il99ll) . An investigation of the basic categories of effective poten- 
tials is presented in Sect. 6 below. 



6 EFFECTIVE POTENTIALS 

The locations of the stationary (stable and unstable) Lagrangian 
points strongly affect the dynamical properties of a rotating galaxy. 
The Lagrangian points and their stability structure the phase space 



at the corotation region. Given a mass distribution, the positions of 
the Lagrangian points are determined by the value of the pattern 
speed Q. p . In the rotating frame we have the conservation of the 
lacobi integral 

Ej = + <l> e ff, (21) 

where <b e /f is the effective potential 

®eff(x,y) = O r (x,y) - ^(x 2 + y 2 ). (22) 

The stationary Lagrangian points at the equilibrium positions are 
defined by 

IT- 1 = I = °- (23) 

ox oy 

This is the case of course, where the system rotates with a sin- 
gle pattern speed. In systems with multiple pattern speeds, multiple 
corotation radii will be defined on the disc. If bar and spirals rotate 
with different Q. p , each component will have its own corotation ra- 
dius. In what follows we take the single pattern speed case as our 
initial assumption. Its validity and its limitation will be examined 
by means of models in subsequent papers. 

Fig. [TO] shows in a compact way how the number, the location 
and the stability of the Lagrangian points change on the disc for the 
various effective potentials. In Fig.llOb we give the positions (x,y) 
of the L\, Li Lagrangian points as a function of Cl p for model A. 
The calculation of the stability of the equilibrium points is done by 
following the sta ndard procedure that can b e found in relevant text- 
books (see e.g. iBinnev & Tremainell2008l , pp. 179-183). The red 
and blue points correspond to the unstable and stable Lagrangian 
points, respectively. At the basis of the cubes in all panels it is plot- 
ted the surface density of the galaxy together with the projections 
of the corresponding Lagrangian points. The projected orange and 
light blue points correspond to the unstable and stable Lagrangian 
points, respectively. Figure [Tob is similar to FigllOb and shows the 
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positions of L 4 , L 5 Lagrangian points in model A. We observe that 
there are ranges of Ci p values for which there are more than one 
Li,Lq,Ls points. Figures [Tob.d andllOfc.f are similar to Fig, UOh.b 
but for the model B and C, respectively. We see that the distribu- 
tion of Lagrangian points have a very similar behavior in models B 
and C. Moreover, the variability of the positions of the Lagrangian 
points with Ci p seems smoother in models B and C. However, even 
in these cases we have ranges of Q p values corresponding to multi- 
ple Lagrangian points (see Fig. II lb). 

Figure QT] presents some characteristic cases of effective po- 
tentials with multiple Lagrangian points. Stable (unstable) La- 
grangian points are plotted with yellow (green) big dots. Fig. [TTb 
shows the isocontours of the effective (total) potential, in model A, 
for Q p = 23.7km sec~'kpc~'. In this figure we observe 5 L\ La- 
grangian points (3 unstable and 2 stable) at the region close to the 
end of the bar. Figure [TTb corresponds to the same model A but 
for Q. p = 21.2km sec~'kpc~'. We observe 3 L5 Lagrangian points 
(2 stable, 1 unstable). Finally, Fig. II It corresponds to model C for 
Q. p = 25.75km sec~'kpc~' and shows 3 L 2 Lagrangian points (2 
unstable, 1 stable). 



7 DISCUSSION AND CONCLUSIONS 

In this paper, we propose three different general models for the po- 
tential of the grand design barred-spiral galaxy NGC 1300, based 
on observations in the near-infrared. The ultimate goal of our work 
is the investigation of the stellar dynamics of this galaxy by means 
of response models and orbital theory. More precisely we want to 
reveal the dynamical mechanisms that shape the observed structure 
of the galaxy. This is done in the forthcoming papers of this series. 

Observations in the near-infrared provide the most reliable po- 
tential models of disc galaxies, since they are excellent indicators 
of the distribution of mass associated to stellar matter. The needed 
additional assumptions refer mainly to the distribution of luminous 
matter in the third dimension (outside the galactic plane) and to the 
distribution of the non-luminous matter (dark halo and central mass 
concentration). 

The three proposed models are: The pure 2D case (Model A), 
where all luminous matter is considered on the galactic plane; the 
thick disc case (Model B) where all luminous matter is considered 
in a 3D disc with a constant scale height; and the "spheroidal case" 
(Model C) in which the major part of the bar is a spherical modified 
Hubble profile. The remaining luminous matter is attributed to a 
3D disc as in Model B. In all these cases the calculation of the 
potential is performed by means of an FFT tra nsform technique, 
that has been introduced in various versions by (Hohl & Hocknevl 
1 19691 ; lOuillen et alii 19941 ; lLaurikainen & Saloll2002l) . 

These three general models can be considered as limiting 
cases in the parameter space of models whose responses and or- 
bital content will be studied. Let us consider an unknown parame- 
ter, e.g. the disc thickness. If we find that we obtain better feedback 
in Model A than in Model B, we will conclude that we have more 
accurate dynamics closer to the thin disc limit. Only the study of 
the dynamics in a large number of effective potentials will show 
which configuration compares best with the observed morphology. 

As we have seen in Section 5, all models are strongly non- 
linear. The relative force perturbation (regarding the axisymmetric 
background) in all models is high reaching to 45%-50% in Models 
B and C and to even higher value 80% in Model A. We note that for 
normal (non-barred) spiral galaxies the models that reproduce suc- 
cessfully the morphology of this type of galaxies have maximum 



force perturbations typically of the order of 5-10% (Patsis et al 



199ll;lLaurikainen & Salol2002; Lauri kainen et alj2004IVoro bvov 



2006) and they are already characterized by nonlinear effects. Non- 
linear effect s are even stronger in larger pertu r bation characteristic 
of bars (see iKaufmann & C ontopoulos 19 9(]; Laurikainen & Said 
l2002l ; lLaurikainen et alj2004l ; lBlock et al.l2004 and others). Thus, 
it is expected that all models will present a strongly nonlinear be- 
havior with extended chaotic regions. Note that we certified the ro- 
bustness of the force perturbation in each model to small changes 
of the various required parameter values (e.g. thickness of the disc). 

The nonlinearity of a model, expressed as the relative force 
perturbation, is an important parameter for studying its dynamics, 
since it refers to the total distribution of matter (luminous and dark) 
in the galaxy. However, the estimation of the potential of a disc 
galaxy alone does not determine its dynamics. The main, and most 
important, parameter is the pattern speed and this cannot be unam- 
biguously determined from observations. The Coriolis forces to- 
tally alter the landscape of the gravitational field, as we can see 
in the various effective potentials. A potential by itself can tell us 
whether or not nonlinear phenomena can play an important role in 
the dynamics of the system. However, for the same potential we 
can have totally different dynamics. By varying the pattern speed 
the equilibrium (Lagrangian) points are shifted and may alter their 
stability character. In our analysis we find cases with multiple La- 
grangian points. This indicates that in real galaxies the existence 
of multiple Lagrangian points can be a common phenomenon. Our 
models give the opportunity to study the dynamics of such realistic 
systems, something that has not been done yet. The effect of the 
variation of the pattern speed on the dynamics of NGC 1300 will 
be presented in the forthcoming papers of this series. 
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